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ABSTRACT 

Distributions of nonthermal surface brightness of supernova remnants (SNRs) con- 
tain important information about the properties of magnetic field and acceleration 
of charged particles. In the present paper, the synchrotron radio, X-ray, and inverse- 
Compton (IC) 7-ray maps of adiabatic SNRs in uniform interstellar medium and 
interstellar magnetic field are modeled and their morphology is analyzed, with par- 
ticular emphasis to comparison of azimuthal and radial variations of brightness in 
radio, X-rays, and 7-rays. Approximate analytical formulae for the azimuthal and ra- 
dial profiles of the synchrotron radio and X-ray as well as the IC 7-ray brightness 
are derived. They reveal the main factors which influence the pattern of the surface 
brightness distribution due to leptonic emission processes in shells of SNRs and can 
account for some non-linear effects of acceleration if necessary. These approximations 
provide observers and theorists with a set of simple diagnostic tools for quick analysis 
of the non-thermal maps of SNRs. 

Key words: ISM: supernova remnants - shock waves - ISM: cosmic rays - radiation 
mechanisms: non-thermal - acceleration of particles 
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\D . 1 INTRODUCTION 

• " Non-thermal images of SNRs are rich sources of information 
^ B about the properties of interstellar magnetic field (ISMF) 
1 behavior and kinetics of charged particles in vicinity of the 
I I strong non-relativistic shocks. Despite of their importance, 

• • . images of SNRs - in contrast to broad-band spectra - are 

^ ' not well studied. 

Synchrotron X-ray brightness profiles were used as di- 
$-H ' agnostic tools for the es timate of the post-shoc k magnetic 
5^ \ field in some SNRs (e.g. iBerezhko fe Volk|[20ul ). Radio az- 
imuthal profiles were used for determ ination of some proper- 
ties of SN 1006 (jPetruk et al.ll2009d . hereafter Paper I) and 
X-ray radial profiles were used to detect the s hock precursor 
in SN 1006 and prove particle acceleration (|Morlino et al.1 
l2010h . 

A detailed approach to modeling the synchrotron im- 
ages of adiabatic SNRs in uniform ISMF and unif o rm in - 
terstellar medium (ISM) is d e veloped by iRevnoldd (I1998T) . 
iFulbright fe Reynolds! (|l990h : iRevnolddl (|l998l . l2004l ~se 
modeled synchrotron maps of SNRs to put constraints on 
properties of accelerated particles. Properties of the inverse- 
Compton (IC) 7-ray maps are invest igated and compared to 
radio images in lPetruk et al.l (|2009al hereafter Paper II). 

The influence of nonuniform ISM and/or nonuniform 



ISMF on the thermal X-ray mor phology of adiabatic SNRs 
are studied in [Hnatvk fe Petrukl ([l999), on the radio maps 
in lOrlando et al.l (120071) and on the syn chrotron X-ray and 
IC 7-ray images in Orlando et al ] (|2010l ). 

All studies of SNR maps assume classic MHD and test- 
particle theory of acceleration. Though they neglect effects 
of the back-reaction of the efficiently accelerated particles, 
they are able to explain general properties of the distribution 
of the surface brightness in radio, X-rays and 7-rays. This is 
because the classic theory, in contrast to the non-linear one, 
is able to deal with oblique shocks, that is vital for synthesis 
of SNR images. 

At present, the theory which considers effects of ac- 
celerated particles on the shock and on acceleration it- 
self is developed for the initially quasi-parallel shocks only. 
One may therefore model the only radial profiles of bright- 
ness, in the rather narrow region close to the shock (in 
order to be certain that obliquity does not introduce 
prominent modifications). Effects of non- linear accelera- 
tio n on the radial profiles of bri g htness are considered 
in Ellison &; Cassam-Chenail (l2005h : ICassam-Chenai' et al.l 



(|2Q05h : iLee et al.1 (|2008l ) : IZirakashvili fe Aharonianl (|2010h . 

Future studies on SNR morphology should take into 
account the NLA effects. Nevertheless, the classic approach 
is able to reveal the general properties of SNR maps de- 
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termined by MF behavior and particle acceleration. Beside 
that, it is important to know the properties of the 'classic' 
images because they create the reference base for investiga- 
tion of the efficiency of NLA effects in the surface brightness 
distribution of SNRs. 

In Paper I, we introduced a method to derive an aspect 
angle of ISMF from the radio brightness of SNR. In Paper II, 
we synthesized radio and IC 7-ray maps and concluded that 
coincidence of the position of the 7-ray and radio limbs is 
not a common case in theoretical models, because different 
parameters are dominant in determintion of the radio and 
7-ray brightness variations . On the other hand, radio, X- 
ray and 7-ray observatio ns ([Petruk e t al. 2QQ9c1; [Miceli et all 
2009; lAcero at al1l2010h show that radio, X-ray and 7-ray 
limbs coincide in SN 1006. As discussed in Paper II, such co- 
incidence might be attributed to a combination of obliquity 
dependences of magnetic field and properties of emitting 
particles, as well as orientation versus the observer. 

In the present paper, we make a step forward extend- 
ing our model of non-thermal leptonic emission of Sedov 
SNRs in uniform medium to the X-ray band. We also derive 
brightness profiles for representative parameters which are 
suited for the comparison with adiabatic SNRs. Moreover, 
we derive analytical approximations of the azimuthal and ra- 
dial profiles of radio, X-ray and 7-ray brightness which are 
extremely useful to demonstrate their dependence on the 
acceleration parameters, magnetic field orientation and the 
viewing geometry; they can also be directly and very easily 
fitted to SNR images to derive estimations for the best-fit 
quantities. While the analytical approximation cannot sub- 
stitute the more accurate numerical simulations, we show 
that they retain enough accuracy to represent an effective 
diagnostic tool for the study of non-thermal SNR shells. 



2 RIGOROUS TREATMENT OF 

SYNCHROTRON AND IC EMISSION OF 
SEDOV SNRS 

2.1 General considerations 

The present paper is continuation of the study in Paper II. 
The model and numerical realisation are therefore similar 
to those used in the Paper II; the reader is referred to this 
paper for more details. In short, in order to investigate the 
prop erties of the leptonic emission of Sedov SNR, we use 
the 
the 



Sedovl Jl 959) solutions for dynamics of fluid as well as 
Reynolds! (| 19981 ) description of the MF behavior down- 



stream of the shock. The use of analytical solutions allows 
us to reduce the computational time considerably. We follow 
IRevnoldsl (| 19981 ) also in calculation of the evolution of the 
distribution function N(E) of relativistic electrons down- 
stream of the shock (see Appendix [Al for details). 

The ISMF orientation versus observer is described by 
the aspect angle O , an angle between ISMF and the line 
of sight. Let us define also the obliquity angle G as the 
angle between ISMF and the normal to the shock, and the 
azimuth angle ip in the projection plane which is measured 
from the direction of the component of ISMF in the plane 
of the sky (see Fig. 1 in Paper II). 

The compression factor for ISMF at the shock front <tb 
is modulated from unity at parallel shock to a = (7+1)7(7 — 



1) at perpendic ular one (where 7 is the adiabatic index), in 
agreement with IRevnoldsl (| 19981 ) . 

At the shock, the energy spectrum of relativistic elec- 
trons is taken to be N(E) = KE~ S exp (-(£/£ max ) a ), 
where E max is the maximum energy of electrons, s the spec- 
tral index, K the normalization, free parameter a regulates 
the rate of the spectrum decrease around the high-energy 
encQ. Evolution of the 'electron injection ability' of the shock 
is represented as K oc V~ b where V is the shock velocity, 
and b is a constant. The variation of the distribution N(E) 
over the surface of the SNR and its evol ution downstream 
of the sho ck are calc u lated as described by IRevnoldsl (|l998h . 
Following IRevnoldsl ([1998), we synthesize the synchrotron 
and IC emission, considering each of the three cases of vari- 
ation of electron injection efficiency with the shock obliquity 
(quasi-perpendicular, isotropic, and quasi-parallel particle 
injection) and each of the three alternatives for time and 
spatial dependence of E m &x (time- limited, loss- limited and 
escape-limited with a gyrofactor 77 as a parameter; it is a 
ratio of the mean free path of a particle to its gyroradius). 
Both the synchrotron and the IC losses are included as chan- 
nels for the radiative losses of relativistic electrons, though 
the IC losses are quite small comparing to synchrotron due 
to large MF considered. 

The surface brightness is calculated by integrating emis- 
sivities along the line of sight within the SNR. The emissivity 
of electrons is given by 



f 

Jo 



N(E)p(E,s)dE 



(1) 



where e is the photon energy and p is the spectral distribu- 
tion of synchrotron or IC radiation power of electron. We 
assume that information about orientation of B inside SNR 
is lost because of turbulence, in practice, we use an average 
aspect angle downstream. The distribution p is calculated 
with the use of the analytical approximation developed by 
IPetrukl (|2009l . see also Paper II). 



2.2 Images 

The resulting synchrotron radio images and the IC 7-ray im- 
ages synthesized by our model have been already presented 
in Paper I and Paper II. Therefore, we present here only the 
X-ray images (see Fig. [1]), adding appropriate references to 
the previous work to let the reader do a quick comparison. 

The pattern of synchrotron X-ray brightness of SNR is 
in general similar to the radio one. In most cases, the bright 
X-ray limbs or other features are located in the same az- 
imuth as in the radio images. The only differences appear 



1 A number of o bservations (lEllison et al.l l200l l200ll ; 

lLazendic et al.l 120041 ; lUchivama et al.l 120031 ) reveal evidence 
of broadening of the high-energy cut-off, i.e. a < 1. Such broad- 
ening should be attributed to the physics of acceleration (|Petrukl 
2006), rather than to the effect of superposition of spectra in 
different co nditions along t he line of sight as it was suggested 
in itially bv | Reynoldsl dl996h . From other side, theoretical model 
of Zirakashvili & Aha ronianl (120071 ) demonstrates that it should 
be a = 2 for the loss-limited models, whereas in practice, taking 
the time evolution into acc ount it could be between 1 and 2 
(|Schure et al.l 120101 ). Also iKang , 
a > 1. 



Rvul (|2010T ) suggest that 
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Figure 1. Maps of synchrotron X-ray surface brightness synthesized for different aspect angles (labeled in the upper left corner of each 
panel). The figure shows the quasi-perpendicular (top), and quasi-parallel (bottom) particle injection models. The model of E mayi is 
time-limited with 77 = 1. Calculations are done for i£ max || = 26TeV, B Q = 30 /iG, s = 2, a = 0.5, b = —3/2. The ambient magnetic field 
is along the horizontal axis. 



due to radiative losses which modify downstream distribu- 
tion of the electrons emitting in X-rays (thus the features of 
brightness are radially thinner) and due to variation of E max 
over the SNR surface. In the radio (see Fig. 4 in Paper II) 
as also in the X-ray band, the remnant shows two symmet- 
ric bright lobes (for <j) — 90°) in all the injection models 
with the maxima in surface brightness coincident in the two 
bands. The maxima are located at perpendicular shocks in 
the quasi-perpendicular and isotropic models (i.e. where B 
is higher), and at parallel shocks in the quasi-parallel model 
(i.e. where emitting electrons are only presented). The lobes 
are much radially thinner in X-rays than in radio because of 
the large radiative losses at the highest energies that make 
the X-ray emission dominated by radii closest to the shock. 

The X-ray morphology of SNR is different for differ- 
ent aspect angles (Fig. [T] cf. with Figs. 4, 5 and 6 of Pa- 
per II for radio and 7-ray images, respectively). In the case 
of quasi-perpendicular injection, the morphology is bilateral 
(two lobes) for large aspect angles (<f> > 60°, i.e. the com- 
ponent of ISMF in the plane of the sky is larger than that 
along the line of sight) and almost ring- like for low aspect 
angles (<j> < 30°; see Fig. [J) with intermediate morphol- 
ogy between 30° and 60°. In the case of quasi-parallel injec- 
tion, the remnant morphology in the radio band is known 
to be bilateral for large aspect angle s and characterized by 
one or two eyes for low asp ect angles ([Fulbright &; Reynolds! 
Il99d : lOrlando et al.1 120071 ). O n the other hand, it is worth 
noting that the remnant morphology in X-rays is in gen- 
eral bilateral for aspect angles (j) > 30° and centrally bright 
for very low angles, indeed a rather limited set of possi- 
ble cases (lower panels in Fig. [J). This happens because the 



non-thermal X-ray emission originates from a very thin shell 
behind the shock making the effect of limb brightening in 
X-rays more important than in the radio band. In addition, 
we note that centrally bright X-ray (and radio) SNRs are 
expected to be much fainter than bilateral SNRs (see lower 
panels in Fig. [TJ and consequently much more difficult to be 
observed. The above considerations may affect the statisti- 
cal arguments generally invoked against the quasi-parallel 
injection (i.e. the fact that th is model produces morphol- 
ogy which is not ob served; e.g. iFulbright &; Reyno lds 1990; 
lOrlando et ai1l2007h . 

Images on Fig. [T] are calculated assuming a dependency 
of -Emax on the obliquity angle which correspon ds to the 
time- l imited model with 77 = 1 as introduced by I Reynolds! 
(|l998h , namely £ max oc ^ m ax(6o)5 , where £ ma x(Bo) is a 
function describing smooth variations of ^ max versus obliq- 
uity, and B Q is the pre-shock ISMF strength (in suc h a model 
no d ependency on the shock velocity is present; Reynolds] 
1998). With this particular choice, ^ max is quite similar 
for different obliquities, namely ^ ma x±/^max|| = 1-3. Larger 
values of 77 always provide ^ m ax±/^ m ax|| > 1 in this model, 
thus the character of azimuthal variation of brightness would 
be similar. 

The efficiency of variation of i^max with obliquity in 
modification of the azimuthal distribution of X-ray syn- 
chrotron brightness depends obviously on the photon en- 
ergy: if the maximum contribution to the emission at given 
photon energy is from electrons with energy much less than 
E max then this effect is negligible. It is useful to introduce 
the reduced photon energy, as e = is/is c (E max \\, B Q ) where 
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Figure 4. Radio (dashed lines) and X-ray (solid lines) radial 
Figure 2. Radial profiles of X-ray surface brightness for different profiles of surface brightness for the adiabatic index 7 = 5/3, 4/3 

photon energies e (marked near respective lines). Calculations are and 1.1. Calculations are done for <f> Q = 0, (p = 0, 6 = 0, efy =3, 

done for O = 0, ip = 0, b = -3/2, e f y = 1. e = 0.3. 
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Figure 3. Effect of radiative losses (represented by efy) on the 
radial profiles of X-ray surface brightness, at two azimuth. The 
reduced fiducial energy e f || = 0.3,1,3 (from below); <f> Q = 7r/2, 
e = 1, 6 = 0. 

v c (E,B) — a (s\ri(f))E 2 B is the synchrotron characteristic 
frequency, c\ — 6.26 x 10 18 cgs, or 

where £kev is the photon energy in keV. If e — 0.29 then 
most of the contribution to the synchrotron X-ray emission 
is from electrons with energy Hi — ^max Figs, n is calcu- 
lated for e = 2.8, i.e. images shown are mainly due to emis- 
sion of electrons with energy few times higher than E max \\; 
the role of variation of E max is therefore clearly visible in 
the images. 

2 For reference: the maximum contribution to synchrotron X-ray 
emission at 3 keV in MF 30 fiG is from electrons with energies 72 
TeV; the maximum contribution to IC 7-ray emission at 1 TeV 
originates from electrons with energies 17 TeV. 



2.3 Brigthness profiles 

In the present subsection, E ma ^ is assumed to be constant 
in time and the same for any obliquity; in addition, isotropic 
injection, s = 2 and a = 1 in the energy spectrum of rela- 
tivists electrons are assumed. 

The radial thickness of features in the X-ray images is 
sensitive to the photon energy: the larger the energy the 
thinner the limbs (Fig. [2}. This is because radiative losses 
E of electrons with energy E is efficient for more energetic 
electrons, E oc E 2 . If e > 0.29 then most of the contribution 
to the synchrotron X-ray emission is from electrons with 
energies E > Em ax where the radiative losses are of the 
main importance. 

An important factor for emission of highly energetic 
electrons is the fiducial energy, which reflects the importance 
of radiative losses in modification of the e lectron distribu - 
tion. It is defined as ef = 637 (B 2 tE max ) 1 (|Revnoldslll998l ). 
or 




Radiative losses are important for ef < 1 and minor for 
ef > 1. Fig. [3] demonstrates how the value of ef affects the 
radial profiles of X-ray brightness: the smaller ef the thinner 
the rim. 

Our model does not include consistently the effects on 
shock dynamics due to back-reaction of accelerated CRs. 
However, we may approach the effect of shock modification 
by considering different values of the adiabatic index 7 which 
is expected to drop from the value of an ideal monoatomic 
gas. In particular, Fig. H] considers the cases of 7 = 5/3 
(for an ideal monatomic gas), 7 = 4/3 (for a gas dominated 
by relativistic particles), and 7 = 1.1 (for large energy drain 
from the shock region due to the escape of high energy CRs) . 
The shock modification results in more compressible plasma 
and, therefore, in the radially-thinner features of the non- 
thermal images of SNRs. A small dist ance between the for- 
ward shock and contact discontinuity ([Cassam-Chenai' et al.1 
I2Q08I : iMiceli et al.ll2009 1 could also be attributed to 7 < 5/3. 
Effect of the index 7 on the profiles of hydrodynamical pa- 
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rameter s downstrea m of the adiabatic shock is widely stud- 
ied (e.g. Sedov 1959): smaller 7 makes the shock compression 
factor higher, 

( 4 for 7 = 5/3, 
a = I 7 for 7 = 4/3, (4) 
[ 21 for 7 = 1.1, 

and the gradient of density downstream stronger (e.g. Ap- 
pendix [Bj, 

( 12 for 7 = 5/3, 
n(f)^f Knr , ftnr=< 25 for 7 = 4/3, (5) 
[ 88 for 7 = 1.1, 

(where n — n/n a , n s is the density immediately postshock 
and r = r/R < 1). In addition to that, the X-ray (and 
also TeV 7-ray) brightness is modified by increased radiative 
losses of emitting electrons. Really, the larger compression 
leads to the higher post-shock MF and thus to increased 
losses, E oc B 2 , which results in turn in the thinner radial 
profiles of brightness. 

It is unknown how the injection efficiency (the fraction 
of nonthermal particles) depends on the properties of the 
shock. We parameterized its evolution as K s oc V~ b where b 
is a constant. Effect of the parameter b on the radial profiles 
of the surface brightness is demonstrated on Fig. \5\ The 
smaller b the thicker the profiles, because there are more 
emitting electrons in deeper layers, which were injected at 
previous times. This property affects the nonthermal emis- 
sion in all bands. However, the effect is less prominent in X- 
rays (and in TeV 7-rays) if radiative losses are quite effective 
to dominate it (see Fig.0 lines for diffe rent photon en e rgies) . 
This is in agre e ment with finding of IParizot et al.l {2006) 
and lVink et al.l (|2006h who showed that, for E > E max , the 
width of the synchrotron limbs does not depend on the shock 
velocity in the loss limited case. Instead, profiles of the radio 
brightness may be used to put limitations on the value of b. 

In a similar fashion, the X-ray and 7-ray radial profiles 
are affected also by the time evolution of the maximum en- 
ergy, i^max oc V q . However, it seems impossible to determine 
q from such profiles because contribution of other factors is 
often dominant. 

An interesting feature of the synchrotron images of 
SNRs is apparent from Fig. [6] The maxima of the radial 
profiles of brightness for different azimuth are located al- 
most at the same distance p from the center of projection, 
for radio and X-rays. Thus, the best way to analyze the 
azimuthal profiles of the surface brightness is to find the po- 
sition p of the maximum for one azimuth and then to trace 
the azimuthal profile of brightness S p (tp) for fixed p. 



3 EXPLICIT APPROXIMATE ANALYTICAL 
FORMULAE FOR SURFACE BRIGHTNESS 
PROFILES OF NON-THERMAL SNR 
SHELLS 

The rigorous model discussed in the previous section is able 
to predict the non-thermal emission of Sedov SNR under 
a variety of conditions. However, in practical applications, 
it may be rather time-consuming to perform an extensive 
parameter space exploration, and the crucial dependence on 



the relevant parameters of the acceleration processes may 
be hidden. 

In order to understand how the properties of MF, elec- 
tron injection and acceleration influence the brightness dis- 
tribution, we derived analytic approximate formulae for the 
azimuthal and radial profiles of the surface brightness of 
adiabatic SNR. In this way, we can easily see what are the 
main factors which determine the pattern of the nonthermal 
images of SNRs, and which of them are mostly responsible 
for the azimuthal variation of the surface brightness and 
which for the radial one. This turns to be extremely use- 
ful in guiding the comparison with real observations. The 
analytical formulae are valid close to the shock only, but 
are adequate to describe azimuthal and radial variations of 
brightness around maxima which are located close to the 
edge of SNR shells. 

3.1 Radio profiles 

Let the evolution and obliquity variation of the electron 
injection efficiency be denoted as V(t)~ b ^(0 o ) and of the 
obliquity variation of MF compression/amplification as 
cb(Oo); for the sake of generality we assume c(B ) and 
cb(Oo) to be some arbitrary smooth functions. Properties 
of the azimuthal and radial profiles of the radio brightness 
is determined mostly by (Appendix IE)) 

Sr((p,Q) OC c(6 , e ff <j>o) ) 0"B (Oo,eff <j>o) ) ^^^Ir (q) (6) 

where 

/ _ 1 1 g (7) 

r VT^¥ + !) ' 

q — p/R, cr is the shock compression ratio, 

Kr = 3b/2 + (2 + s)n ad + l/cr + s, (8) 

ftad is close to unity for 7 =1.1 -=-5/3 (Appendix [Bj . The 
effective obliquity angle G , e ff is related to azimuth ip and 
aspect <j)o as 

cos9 ,eff (<p, 4>o) = cos sin O , (9) 

the azimuth angle is measured from the direction of ISMF 
in the plane of the sky. Eq. ([6} is a generalisation of the 
approximate formula derived in Paper I. 

Eq. ((6} shows that the azimuthal variation of the radio 
surface brightness S e (ip) at a fixed radius g of projection, 
is mostly determined by the variations of the magnetic field 
compression (and amplification, if any) <tb and by the vari- 
ation of the electron injection efficiency q. The radial profile 
Sip(g) oc I t (q) is determined mostly by a, b and s. Adia- 
batic index 7 affects the radial and azimuthal profiles mostly 
through the compression factor a = (7+1)/ (7 — 1) because 
ftad weakly depends on 7. 

3.2 Synchrotron X-ray profiles 

Let us assume that the maximum energy is expressed as 
^max(Bo,t) oc V (t) q S m a^ (Bo) , where q is a constant and 
£max(B ), for the sake of generality, is some arbitrary func- 
tion describing the smooth variation of E m &x versus obliq- 
uity. The synchrotron X-ray brightness close to the forward 
shock is approximately (Appendix |C|) 
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Figure 5. Evolution of injection efficiency (represented by b) and radial profiles of surface brightness in radio (left), X-rays (middle) 
and 7-rays (right). The parameter in relation K s oc V~ h is b = —3/2, 0, 2 (from above). Other parameters are cf) = 7r/2, <p = (for radio 
and X-rays) and cp = tt/2 (for 7-rays), e-f n = 1. X-ray and 7-ray profiles are shown for two photon energies: where the most contribution 
is from electrons with E = £ max (dashed lines) and three times larger (solid lines); in case of X-rays, solid lines correspond to e = 1. 




Figure 6. Radial profiles of the radio (left), X-ray (middle) and 7-ray (right) surface brightness of SNR for different azimuth: 
0, 7r/6, 7r/3, 7r/2 (marked near lines). The aspect angles are tt/2 (solid lines) and tt/3 (dashed lines); b = —3/2, e f || = 1, e = 1, 7- 
ray profiles are for photons with energy 0.l£ 7 (£^ max )- 



5 X (<£, Q) ^(^)ctb(^) (s+1)/2 exp 
where 

irx = It(q)Ix(<P,Q) 

with 

1 



em(y) 

^-rnax (<^) 



LAv, e)(io) 



(ii) 



Co. 
<-rna 



T(^ r + 2) 



ftr + 1 



1 _ *r(« r + l) Kr _|_ 2 



The parameter 



2e 



2 



(12) 



(13) 



is responsible for the adiabatic (the first term) and radiative 
(the second term) losses of emitting electrons and the time 
evolution of E max on the shock (the third term) . The reduced 
electron energy which gives the maximum contribution to 
emission of photons with energy e is 



(o.29a B ) 



1/2 



(14) 



it varies with obliquity (since MF varies; electrons with dif- 
ferent energies contribute to the synchrotron emission at e) . 
Parameters ip, e m , £max depend on B , e ff and, therefore, on 
the aspect angle <j) and the azimuth angle (p. 

If e m 1 then Eq. (|10p for the X-ray brightness trans- 
forms to Eq. (0 for the radio brightness. 



The thickness of the hard X-ray radial profile is used 
to estimate the post-shock st rength of MF in a number of 
SNRs (e.g. iBerezhko fe Volldl2004h . The absolute value of 
MF is present in Eq. JTDJ through e and e f , Eqs. J2]), © 
which appear in ip and e m , Eqs. (|13[) , (|14|) . In both cases, i3 s 
is in combination with ^ max (thus, the value of the electron 
maximum energy may affect the estimations). The idea of 
the method bases on the increased role of losses in X-rays 
due to larger MF, i.e. on the role of the second term in 
ip, Eq. ([T3|) . Really, the influence of k t (i.e. of s and b) is 
minor in X-rays (middle panel on Fig. [5]), if radiative losses 
affect the electron evolution downstream of the shock (i.e. 
for e >0.29, e f <1). The role of the first and the third 
terms in ip are also minor in most cases (q = for the time- 
limited and escape-limited models and unity for the loss- 
limited one) because the second term > 10. However, the 
adiabatic index makes an important effect on the thickness 
of the profile, mostly through a which appears in as and in 
7 X . Being smaller than 5/3 (that is reasonable especially in 
the case of efficient acceleration, which is actually believed 
to be responsible for the large MF) , the index may compete 
to some extend the role of losses, used in the method for 
estimation of MF (see e.g. Fig. 2} that might lead to smaller 
estimates of MF strength. 
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3.3 IC gamma- ray profiles 

The IC 7-ray brightness may approximately be described as 
(Appendix [D|) 



Sic(<p, g) oc c;((^)exp 
where 



1 



-<r(«; ic + l) 



(15) 



x 



g2 cr(ft ic + 1 
1 



-a(« ic +2) ^ + x 



££ax 1 - + + 2 

K ic = 36/2 + (2 + s)K ad + 1/cr - 1. 



(16) 



(17) 



The expression for ip is the same as ([T3]) but e m is different: 

£ 1/2 

£m = 2(^T)V2 7max|| ' ( 18 ) 

where e is the 7-ray photon energy T the temperature of 
the seed photon field, 7 max the Lorentz factor of electrons 
with energy -c/max 5 or 

- = 1 - 66, K^)" 1 (liv) 1/2 (2iK)" 1/2 - < 19 > 



Eq. (p"5|) is a generalisation of the approximate formula de- 
rived in Paper II. 

The azimuthal variation of the IC 7-ray brightness de- 
pends mostly on the injection efficiency The role of variation 
of i^max is prominent only if obliquity dependence of injec- 
tion is not strong. Parameter a, being smaller than unity 
results in smaller azimuthal contrasts of synchrotron X-ray 
or IC 7-ray brightness comparing to model with purely expo- 
nential cut-off in N(E). The radial distribution of IC bright- 
ness is determined mostly by a, e m , ef , b and, to the smaller 
extend, by s and q. 

3.4 Accuracy of the formulae 

The approximations presented above do not require long and 
complicate numerical simulations but restore all the proper- 
ties of nonthermal images discussed in the previous sections 
as well as in Papers I and II, including dependence on the 
aspect angle. Therefore, they may be used as a simple diag- 
nistic tool for non-thermal maps of SNRs. 

The formulae are rather accurate in description of the 
brightness distribution close to the shock. They do not rep- 
resent centrally-brightened SNRs. Instead, they may be used 
in SNR shells for those azimuth <p where e m <1 and ef > 0.1, 
in the range of g from (1 — 2Ag m ) to 1, where Ag m = 1 — g m , 
g m is the radius (close to the shock) where the maximum 
in the radial profile of brightness happens. Approximations 
are compared with numerical calculations in respective Ap- 
pendices and their applicability is discussed in details on 
example of the IC emission in Appendix ID II 



4 DISCUSSION 

Analysis of azimuthal profiles of brightness in different bands 
allows one to put limitations on models of injection, MF, 



E max . In most cases, the best way to estimate the azimuthal 
variation S p ((p) of the surface brightness is the following. An 
approximate radial profile S^(p) of the brightness should be 
produced for azimuth where the largest losses occur (i.e. 
where ef oc (^maxCs) is smaller; e.g. at if = 90°). This 
allows us to find g m which should be used in S p ((p) in order 
to estimate the azimuthal variation of brightness. g m keeps 
us at maxima in the radial brightness profiles for different 
azimuth (Fig. [6}. 

Energy of electrons evolve downstream of the shock as 
E(p) — EiSadSr&d, where Ei is initial energy, a = a/ 'R, a 
the Lagrangian coordinate. Adiabatic and radiative losses 
of electrons in a given fluid element a are represented by 
functions £ a d(a) < 1, £rad(«) < 1 respectively (Appendix 
[A)) . Close to the shock, they are approximately £ a d ~ a K&d 
(ftad depends on the adiabatic index 7 only and is close to 
unity for 7 = 1.1^5/3), £ rad « a 5(T B e -/( 2e f||) (Appendix©. 
The latter expression clearly shows that the fiducial energy 
€f is important parameter reflecting the ' sensitivi ty' of the 
mode l to the radiative losses, as it is shown by iRevnoldsl 
(1998{): the larger the fiducial energy the smaller the radia- 
tive losses. In fact, £ ra d = 1 means no radiative losses at 
all. Another fact directly visible from this approximation is 
that radiative losses are much more important at the per- 
pendicular shock (where <tb = cr) than at the parallel one 
(where <tb = 1)- In addition, the radiative losses depends 
rather strongly on the index 7: a\ — 16 for 7 = 5/3 but 
al = 49 for 7 = 4/3. 

Note, that in the analysis above, the difference be- 
tween the parallel and perpendicular shocks are only due 
to the compression factor <tb which may be treated as 
"compression-plus-amplification" factor. Therefore, our con- 
sideration may also be applied in case of the very turbulent 
and amplified pre-shock field (when information about obliq- 
uity is lost) once this factor is known. 

Our approximations reflects also the general 'rule' for 
IC emission: there is less IC emission where MF is stronger. 
Namely, the azimuthal variation is Ii c (<p) oc 1 — const x 
ip((p) I 'T^Lp) " with ip oc (Jb(^p) 2 '- emitting electrons disap- 
pears toward the shock with larger ctr because MF strength 
is a reason of higher losses there. Similar dependence on 
<7b is for X-rays, Eq. (|12|) ; it is however dominated by the 
increased term S x ((p) oc ctb(^) 3 ^ 2 - 

TeV 7-ray image of SN 1006 demonstrates good corre- 
latio ns with X-ray imag e smoothed to the H.E.S.S. resolu- 
tion ()Acero at al.ll2010h . We mean here both the location 
and sizes of the bright limbs. Let us consider the polar-caps 
model of SN 1006. The shock is quasi-parallel around the 
limbs, MF azimuthally increases (in w 4 times). The injec- 
tion efficiency decreases (in > 10 -3 times) out of the limbs; 
the number of electrons emitting in X-rays and TeV 7-rays is 
therefore dramatically low at perpendicular shock compar- 
ing to the parallel, that is in agreement with no TeV emission 
at NW and SE regions of SN 1006. However, the azimuthal 
sizes of the limbs in X-rays and 7-rays is expected to be 
different, in the polar-caps model: they should be larger in 
X-rays. Really, azimuthal variation of £ ma x is smaller than 
variations of q and cr 3 / 2 ; therefore, from Eqs. ([TO) and ([T5]h 
azimuthal variations of brightness are mostly S x oc c;cr 3//2 
while S\ c oc q (see Fig. [71 for a comparison of remnant mor- 
phologies in the radio, X-ray and 7-ray bands). We hope that 
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Figure 7. Radio (left), X-ray (center) and 7-ray (right) images of the polar-cap SNR, assuming quasi-parallel injection. Model parameters 
are the same as in Fig. ^ Different azimuthal and radial sizes of limbs in various bands are clearly visible. 



future observations allow us to see if there is a difference in 
azimuthal sizes of the limbs in various bands. 

How back-reaction of accelerated particles may modify 
nonthermal images of SNRs? Our formulae can restore some 
of these effects. In our approximations, s, in general, is al- 
lowed to vary with E, e.g. to be s(E) = s + 5s(E). The index 
s reflects the 'local' slope of the electron spectrum around 
e m . Therefore, if s(E) 7^ const, the index s has to vary with 
azimuth because e m varies, Eq. (|14|) . Generally speaking, 
such approach allows one to estimate the role of the nonlin- 
ear 'concavity' of the electron spectrum in modification of 
the nonthermal images. However, we expect that this effect 
is almost negligible because Ss (E) is very slow function, at 
least within interval of electron energies contributing to im- 
ages. Another effect of efficient acceleration consists in the 
adiabatic index 7 smaller than 5/3. Our approximations are 
written for general 7. Namely, the index 7 affects S((p,g) 
through a. Cosmic rays may also cause the amplification 
of the seed ISMF. In our formulae, ctb(0 o ) represents the 
obliquity variation of the ratio of the downstream MF to 
ISMF strength, B S /B Q . Therefore, it may account for both 
the compression and amplification of ISMF; for such pur- 
pose, <tb should be expressed in a way to be unity at parallel 
shock. 



5 COMPARISON WITH OBSERVATIONS: A 
WORKING EXAMPLE 

The main reason for the derivation of the explicit approx- 
imations of Sect. [3] is to highlight the factors which are 
the most efficient in the formation of the pattern of surface 
brightness. However, sometimes the analytical formulae may 
help in estimates some of the remnant parameters. We would 
like to present two examples. Let us consider SN 1006, s — 2 
(jMiceli et al.ll2009h . 7 = 5/3. 

The injection efficiency is isotropic if one assumes that 
SN 1006 evolves in the uniform ISMF and uniform ISM (Pa- 
per I) . The best-fit value of the aspect angle found from the 
approximate Eq. @ is 68° ± 4.0° (Fig. [8} while the detailed 
numerical calculations give 70° =b 4.2° (Paper I). 

The same approximate formula shows that the radial 
profile of radio brightness depends only on the value of b 



1 




q 2 I i i I i i I i i I i i 

"-90 -45 45 90 

cp, degrees 

Figure 8. Azimuthal profiles of the radio brightness of SN 1006. 
Experimental data and the best-fit theoretical profile calculated 
from detailed model (dashed line) are from Paper I. Green solid 
line represents the best-fit profile obtained from the approximate 
analytical formula (|6|). 

(which shows how the injection efficiency evolve with the 
shock velocity). Unfortunately, the differences between pro- 
files for b = — 1 and b — 1 are comparable with accuracy 
of the approximate formulae and of the experimantal data; 
therefore, the approximation may not be used for estima- 
tions of b. 

The sha rpest radia l profi le of X-ray brightness (from 
Fig. 4A in lLong et al.1 l2003h was used to estimate the 
strength of the post-shock MF in ~ 100 /xG ([Berezhko et al.l 
2003). Fig. [9] shows approximate profiles for three values of 
MF in comparison with the Chandra profile and detailed 
numerical simulations (Petruk et al., in preparation). One 
can also see from the approximation that B ~ 100 /xG is the 
most appropriate value in agreement with the value found 
in the literature. 



6 CONCLUSIONS 

The present paper extends analysis of properties of the 
surface brightness distribution of spherical adiabatic SNRs 
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p/R 



Figure 9. Radial profiles of the hard X-ray brightness of SN 1006. 
Experimental data (from Fig. 4A in lLong et al.ll2003h are shown 
by histogram. Profiles given by approximation (|10|) for B = 
50, 100, 150 /iG are shown by blue, green and red lines respec- 
tively. Other parameters: e = 1.2 keV and ^break|| — 150 eV, 
cf) = 68°, if = 70°, £ max = 2.9, q = 0, b = 0, a = 1. 



started in Paper I (radio band) and Paper II (IC 7-rays) to 
the nonthermal X-rays. It also generalizes the method of ap- 
proximate analytical description of the azimuthal and radial 
profiles of brightness introduced in these papers. 

Synchrotron images of adiabatic SNR in X-rays are 
synthesized for different assumptions about obliquity vari- 
ations of the injection efficiency, MF and maximum energy 
of accelerated electrons. We analyze properties of these im- 
ages. Different models of electron injection (quasi-parallel, 
isotropic and quasi-perpendicular) as well as models of the 
electron maximum energy (time- limited, loss- limited and 
escape- limited) are considered. 

The azimuthal variation of the synchrotron X-ray and 
IC 7-ray brightness is mostly determined by variations of 
(JB and E max , of the radio brightness by <r and <tb only. 
In general, higher B increases X-ray and decreases IC 7-ray 
brightness. Really, higher MF is a reason of larger losses of 
emitting electrons (i.e. decrease of their number) and thus of 
the smaller brightness due to IC process. In contrast, X-rays 
are more efficient there because S x oc B 3 ^ 2 . 

The radial profiles of brightness depend on a number of 
factors. It is quite sensitive to the adiabatic index: 7 < 5/3 
makes plasma more compressible. Therefore, the brightness 
profile is thinner due to larger compression factor, larger 
gradient of density and MF downstream of the shock and 
larger radiative losses. 

The role and importance of various factors on the sur- 
face brightness in radio, synchrotron X-rays and IC 7-rays 
are demonstrated by the approximate analytical formulae. 
They accurately represent numerical simulations close to the 
shock and are able to account for some non-linear effects of 
acceleration if necessary. This makes the approximations a 
powerful tool for quick analysis of the surface brightness 
distribution due to emission of accelerated electrons around 
SNR shells. The application of the approximate formulae to 
the case of SN1006 yields measures of the aspect angle and 
the post-shock MF in good agreement with more accurate 
analysis found in the literature. 
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APPENDIX A: EVOLUTION OF THE 
ELECTRON ENERGY SPECTRUM 
DOWNSTREAM OF THE ADIABATIC SHOCK 

Relativistic electrons evolving downstream of the shock suf- 
fer from adiabatic expansion and radiative losses due to syn- 
chrotron and inverse Compton processes. Electrons are con- 
sidered to be confined to the fl uid element w hich removed 
them from the acceleration site ([Reyno lds 1991) 

Let the fluid element with Lagrangian coordinate a was 
shocked at ti me ti. If energy of electrons was E[ at ti, it 
becomes later (Reynolds 1998) 



fi(a) ~ V3 +I(a)Ei/Ef 



(Al) 



where the first summand in the denominator reflects adi- 
abatic losses, the second one is due to radiative losses, 
n — n/n s , index 's' means 'immediately post-shock', the 
fiducial energy for parallel shock Ef = 637 /(B 2 ^ ;S m t) cgs. 
The effective magnetic field is B 2 ^ — B 2 + E^ M r , B(a) take s 
into account its evolution downstream ([Reynolds! [l998), 
Bcmb =3.27 /nG strength of magnetic field with energy den- 
sity equal to energy density of CMBR. Bcmb is in troduced 
in or der to account for the inverse- Compton losses ([Reynolds! 
1998), therefore it is constant everywhere. The synchrotron 
channel dominates inverse- Compton losses if B s ^> Scmb- 

The dimensionless function X accounts for evolution of 
fluid during time from tj to t; it was initially defined as 
integral over time (Reynolds 1998). In case of Sedov shock, 
X may be written in terms of spatial coordinate that is more 
convenient for simulations than original representation in 
terms of time. Namely, for uniform ISM: 



X(a, Go, d) 



5ai 



3/2 ,) 

x 1 B ( 



2n(o)V3 J a 
Eq. (|A1[) results in relations 



-© 2 ^© 1/3 ^ (A2) 



E { = 



dEi 



dE 



£ad<^rad ^ a( ^^rad 



(A3) 



where the adiabatic and radiative losses are represented by 



-ad — n( a 



,V3 



£ rad = l-X(a,e o )E/E t . 



(A4) 



Shocks of different strength are able to accelerate elec- 
trons to different E max . Let E ma ^ oc V q where V is the 
shock velocity, q = 1,2,0 for loss-li mited, time-lim ited and 
escape- limited models respectively (Reynolds 1998). If shock 
accelerates electrons to ^ max at present time t, then, at some 
previous time ti when fluid element a = R(ti) was shocked, 
the shock was able to accelerate electrons to 



-E'max(ti) — E n 



V v(t) J 



-3q/2 



(A5) 



where a = a/R(t) and we used Sedov solutions. The obliq- 
uity variation of the maximum energy is given by E max = 
(Go) with £max(G ) independent of time. 



^max, || £max 



3 Such an assumption means that we consider the advection 
only, not the diffusion. However, for X-ray emitting electrons, the 
lengthscales for advec tion t, a 

'advAdif — 0.< 



dv an d diffusion Z^if are comparable: 
1.1 (Ballet 2006). Therefore, the results of the 



Let us assume that, at time ti, an electron distribution 
has been produced at the shock 



N(Ei,ti) = ^ s (ti,G )^r s exp 



<(ti,G ) 



(A6) 



where a is constant. The obliquity variation of K s is given by 
K S (G ) — i^ s ||/C(G ) with /C(G ) also independent of time. 
Conservation equation 

N(E,a,t) = N(E i ,a,t i) ^§ (AT) 

(where a = n s /n ) and continuity equation n Q (a)a 2 da = 
n(a,t)r 2 dr shows that downstream 

N(E, a, t) = K(a, t, G G ) E~ s S rad (a, E, G ) s ~ 2 

Ea^ 2 



x exp 



(t) ^ad(a) S rdL d(a,E) ^max(Go) 



(A8) 



with K{a,t) = Ks^fiS*- 1 . If K s oc V~\ then evolution of 
K is self-similar downstream 

K(a) = K(a,t)/K s (t) 

Therefore, in general, 

K(a, t, Go) = K sll (t)JC(Q )K(d 



(A9) 



(A10) 



where the profile K(a) is independent of obliquity. Note, 
that K(a, t, G Q ) is not affected by the radiative losses, there- 
fore it behaves in the same way also for radio emitting elec- 
trons. Once s is close to 2, the radiative losses influence 
the shape of N(E) mostly through the exponential term in 
Eq. (|A8[) . In other words, they are effective only around the 
high- ener g y end of the electron spectrum as it is shown by 
iRevnoldsl (jl998h . 

The above formulae are also valid if the spectral index s 
depends on E, e.g. s(E) = s + 5s(E) like it would be in case 
of the nonlinear acceleration. In addition, no specific value 
of the adiabatic index 7 is assumed here. It influences the 
downstream evolutio n of relativis tic electrons through n(a) 
which depends on 7 (|Sedovlll959h . 



APPENDIX B: APPROXIMATIONS FOR 
DISTRIBUTIONS OF SOME PARAMETERS 
BEHIND THE ADIABATIC SHOCK 

Let us find approximations for dependence of some param- 
eter X = X/X(R) on the Lagrangian coordinate a = a/R 
downstream close to the adiabatic shock. We are interested 
in approximations of the form 



X(a) « a K 
where, by definition, 
a dX* (a) 



X* (a) da 



dhiX*(a) 
d\na 



(Bl) 



(B2) 



and star marks the dependence given by the Sedov solution. 
This approach yields for density 

5 7 +13 



n{a) : 



(7 + 1) 2 ' 



(B3) 



present paper are robust. 



for the relation between Eulerian and Lagrangian coordi- 
nates 
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,1/a 



rtta>~, r, « (l/a)a (1/CT)_1 (B4) 
where the shock compression factor is 

, = I±1. (B5) 

Note that the density distribution in Eulerian coordinates is 
much more sensitive to 7 fTable lBT]) : 

5 7 +13 



n{r) w r Knr , ^ nr = an na = 



(B6) 



(7+l)( 7 -l)- 
Magnetic field is approximately 

5~^ (0 °' eff) , (B7) 
cos 2 Go + /3±cr 2 sin 2 G Q 



/3(6 ,eff 



cos 2 Go + a 2 sin 2 O 
37 + 11 



7 + 1' 



Pi 



(7+1) 



(B8) 
(B9) 



Approximation for normalization K follows from the defini- 
tion (Appendix [A]) 



K = a so/ 'n 



36/2-(2+s)/3 



(BIO) 



and approximation for n. 

Adiabatic losses are accounted with S a d(a) which is de- 
fined by (|A4[) . Its approximation is therefore 



57 + 13 
3(7 + l) 2 



(Bll) 



it is valid for f > 0.8 with error less than few per cent. The 
value of ftad is close to unity for 7 =1.1 + 5/3 (Table IB ip . 

In order to approximate £ ra d defined by (|A4[) , we sub- 
stitute (|A2[) with approximations for n and B. Then we use 
the property 



lim 



cy 



(B12) 



o-n \ f(a) t 

for function of the form f(a) = 1 — ca x (l — a y ). In this way, 

5a!(e o )£ 



£ ra d(a, E 1 ) « a K 



Krad 



2E fl 



(B13) 



This expression is good for r > 0.94, with error of few per 
cent. It depends on 7 through a in <tb which is ( Reynolds! 
I1998T ) 



CTB = 



1 + a 2 tan 2 6 x 1/2 
1 + tan 2 Go 



(B14) 



Note that dependence on the absolute value of the mag- 
netic field strength B s \\ is present in (|B13[) : E^\ oc 5^ 2 - This 
approximation clearly shows that the fiducial energy Ef is 
important parameter reflecting the 'sen s itivity ' to the ra- 
diative losses, as it shown by I Reynolds! (|l998h : the larger 
the fiducial energy the smaller the radiative losses. In fact, 
£rad = 1 means no radiative losses at all, see (|A4[) . An- 
other fact directly visible from Eq. (|B13[) is that radiative 
losses are much more important at the perpendicular shock 
(ctb = cr) than at the parallel one (<tb = 1)- Radiative losses 
depends rather strongly on the index 7: cr B = 16 for 7 = 5/3 
but <t% = 49 for 7 = 4/3. 

The values of parameters in approximations for different 
adiabatic index 7 are presented in Table IB1I 
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Table Bl. Parameters in approximations 



Expression 



7 = 5/3 7 = 4/3 7 = 1.1 



_ 5 7 + 13 

= 57 + 13 

* nr " (7 + 1X7-1) 
7 + 1 

a = 

7-1 
37 + n 

^ad — ^na/3 



3 


3.6 


4.2 


12 


25 


88 


4 


7 


21 


1.5 


1.7 


1.9 


2.2 


2.8 


3.2 


1 


1.2 


1.4 



APPENDIX C: APPROXIMATE FORMULA 
FOR THE AZIMUTHAL VARIATION OF THE 
SYNCHROTRON X-RAY SURFACE 
BRIGHTNESS IN SEDOV SNR 

A formula obtained here may be useful in situations when 
an approximate quantitative estimation for the azimuthal 
variation of the synchrotron X-ray surface brightness is suf- 
ficient. 

1. The emissivity due to synchrotron emission is 



S 



N(E)p(E,s)dE 



(CI) 



Spectral distribution of the synchrotron radiation power of 
electrons with energy E in magnetic field of the strength B 
is 



P(E, 



v/3e 3 Bsin( 



(C2) 



where v is frequency, v c — c\BE 2 the characteristic fre- 
quency. Most of this radiation is in photons with energy 
s p = 0.29/iz^c- In the 'delta- function approximation', the spe- 
cial function F is substituted with 



00 

F (v)= 5 {v c - Q ™)l F{x)dx 



where 



F(x)dx= 

9V3 



With this approximation, (|C1[) becomes 
47re 3 sin 6 e 1/2 B 1/2 



9m e c 2 0.29ci 7 V/2 



N(E n 



(C3) 



(C4) 



(C5) 



where E m is the energy of electrons which give maximum 
contribution to synchrotron emission at photons with energy 
e: E m = £ 1/2 (0.29/iciB)" 1/2 . 

2. Let the energy of relativist ic electrons is E in a given 
fluid element at present time. Their energy was Ei at the 
time this element was shocked. These two energies are re- 
lated as 



E — EiS a dS Tl 



(C6) 



© 2010 RAS, MNRAS 000.HHT51 



12 Petruk 0. et al. 



where £ a d accounts for the adiabatic losses and £ ra d for the 
radiative losses (Appendix [Aj). There are approximations 
valid close to the shock (Appendix [Bj: 



Sri 



-5<T^E/2E tA 



(C7) 



where a = a/R, a is Lagrangian coordinate of the fluid ele- 
ment, £7 fj || is the fiducial energy for parallel shock, ^ a d de- 
pends on 7 and is given by (|B11[) : ^ a d = 1 for 7 = 5/3 
(for other 7 see Table IB1[) . Th e factor ap re presents com- 
pression in the classical MHD (|Revnoldjll99a ) but may be 
interpreted also as amplification-plus-compression factor. In 
the latter case, it should be written in a way to be unity at 
parallel shock. 

The downstream evolution of K in a Sedov SNR is (Ap- 
pendix [A| 



K oc q(e o )K(a) 



(C8) 



where c; is injection efficiency. With the approximations (|C7|) 
and s close to 2, the distribution N(E) may be written from 
[8]) as 



N(E, Bo) oc s(0 o )K(a)E- s exp 



Ea 



H^max(Oo) 



(C9) 



where 



ip(E,e ) = n ad + 



5a B (e o ) 2 E 
2E f I, 



34 
2 



(CIO) 



and s is allowed to vary with E. 

3. Let us consider the azimuthal profile of the syn- 
chrotron X-ray brightness S e at a given radius g from the 
center of the SNR projection. 

Like in Paper II, we consider the 'effective' obliquity 
angle G , e ff which, for a given azimuth, equals to the obliq- 
uity angle for a sector with the same azimuth in the plane of 
the sky (see details in Paper II). The relation between the 
azimuthal angle ip, the obliquity angle o , e ff and the aspect 
angle <p Q is as simple as 



C0s9 ,eff 



cos ip sin c 



(Cll) 



for the azimuth angle ip measured from the direction of ISMF 
in the plane of the sky. 

The surface brightness of SNR projection at distance g 
from the center and at azimuth <p is 



S(q, ip) 



q(a) 



rr^G 



a(g) 



(C12) 



where f a is the derivative of f(a) in respect to a. The az- 
imuthal variation of the synchrotron X-ray brightness is ap- 
proximately 

s x oc ^(e , eff ) ( 7B(e 0)eff )( s + 1 )/ 2 / rx (e , eff , e) 

^msfc ©o,eff) 



x exp 



c(©o,eff ) 



(C13) 



where 



KB (s+l)/2 f 



a(g) 



reflects the dependence on p, E ms is E : 

E ms (£, o ,eff) 



for B = l: 

1/2 



(C15) 



0.29/ici£ o a B (e o ,eff) 

Note, that E ms oc £ 1//2 , i.e. S e depends in our approximation 
on the energy e of observed X-ray photons. 

4. Let us approximate 7 rx . First, we use the approxima- 
tions a ~ f a , KB^ s+1 ^ 2 fa « a Kr /cr, which are valid close 
to the shock (Appendix [Bj , a is the shock compression ra- 
tio. Next, we expand r/^/r 2 — g 2 in powers of the small 
parameter (r — 1) and consider the only first term of the 
decomposition: 

f 1 



(C16) 



The exponential term in the integral expands in powers of 
the small parameter (1 — a): 

exp (- Xl (a~ X2 - 1)) « 1 - xix 2 (l - a). (C17) 

In addition, E ms is used instead of E m . 

Close to the shock, the integral of interest is therefore 

Irx(<P,Q)™Ir(p)I X (<P,Q) (C18) 

where 



£0-(«r + l) 



Kr + 1 



Co, 



-o-(« r + 2) 



+ 1 



T(«r + 1) Kr + 2 



The parameter 



Ip = K a d + 



5cr B e m 
2e f „ 



3q 
2 



(C19) 



.(C20) 



(C21) 



is responsible for the losses of emitting electrons and the 
time evolution of E ma ^ on the shock. The value of ft a d is 
rather close to unity for 7=1. l-j-5/3 ( Table iBl]) ; unless 
radiative losses (the second term in ip) are negligible, one 
may use K a d ~ 1 for any 7. Other parameters are 

E ( e \ 1 / 2 

£m = E^jt = (o29^j ' (C22) 



is given by Eq. (|B8[) . 

36 2 + s s + L 1 



(C23) 



(C14) 



x exp 



( a -a^m)s-a/2_^ 



da, 



Parameters ip, e m , ctb, £ ma x and /3 depend on G , e ff and 
therefore on the aspect angle <p Q and the azimuth angle ip. 

The parameter ft reflects differences between MF dis- 
tribution downstream the shock of the different obliquity. It 
varies from /3\\ at parallel shock to /3± at perpendicular one, 
Eq. (|B8[) . In the approximate formulae, it appears in the 
combination /3/2; the role of f3 G is minor in mod- 

ification of the approximate azimuthal and radial profiles. 
Therefore, in order to simplify the approximation, we may 
take f3/2 « 1. 

The index s in (|C9|) . in general, is allowed to vary with 
E, e.g. to be s(E) = s + Ss(E). In our approximation, due 
to (|C3[) , s reflects the 'local' slope of the electron spectrum 
appropriate to e m . Therefore, if one assumes s(E') 7^ const, 
the index s(e m ) may vary with azimuth because 6m varies, 
Eq. (IC221) . 
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5. The final formula is 



exp 



V> = >W + 

g; e f ||)(C24) 



5cr|e n 
2e f ,i 



3g 
2 ' 



r l/2 



depends on £> and 6f|| 
C24l) 



where only 7 : 

The formula Eq. (|C24|) gives us the possibility to ap- 
proximate both the azimuthal and the radial brightness pro- 
files of X-ray brightness for g close to unity. It may be used 
(with a bit larger errors compared to the case of IC emis- 
sion; Fig. lCU cf. Fig. lDl[) , for those azimuth ip where e m < 1 
and ef >0.1, in the range of g from 1 — 2Ag m to 1, where 
Ag m — 1 — g ni . g m is the radius where the maximum in 
the radial profile of brightness happens. We have in mind 
the maximum which is close to the shock, say g m > 0.95; 
therefore, in order to determine £ m , one should look for the 
azimuth with the largest radiative losses. This is discussed 
in details on example of the IC emission in Sect. iDl 

Adiabatic index 7 affects the approximation through cr, 



Cm = 



2(£T)V2 7n 



3b 2 + s 



The final formula is 



Sic(<p, g) oc q((p) exp 



<(*>) 



Iic(<f,Q] Cf||) 



(D5) 



(D6) 



(D7) 



(D8) 



It gives us the possibility to approximate both the azimuthal 
and the radial brightness profiles for g close to unity. 



Dl Accuracy of the approximation 



APPENDIX D: APPROXIMATE FORMULA 
FOR THE IC GAMMA-RAY SURFACE 
BRIGHTNESS IN SEDOV SNR 

In Paper II, we have developed an analytic approximation 
for the azimuthal variation of the surface brightness of Se- 
dov SNR in 7-rays due to the inverse- Compt on process, 
for regions close to the forward shock. The approximation, 
Eq. (11), in the cited paper accounts to zeroth order. How- 
ever, like in the case of the X-ray brightness, the fall of the 
7-ray emissivity downstream of the shock is quite strong in 
case of the efficient radiative losses of electrons. Therefore, 
in such cases of the efficient losses we need to consider the 
next order of approximation. 

Adopting the approach from the Appendix \C\ to IC 
emission (see also some details in Paper II), we come to 
the approximation 



Sic(<p) oc c(0 o ,eff) exp 



E m (e) 



|| ^max(0o,eff ) 



/ic(e. 



where the energy Em of electrons which gives maximum con- 
tribution to IC emission at photons with energy s is (e.g 
IPetnildB) 



m e cV/ 2 
m 2(fcT)V2 ' 

T is the temperature of the black-body photons. 
The factor 



(D2) 



h 



Krr&da 



/~2 



is approximately 



Q 2 



exp 



(a 



-aip 



1 



g^ 



c + l) 



Col 



Kic + 1 



prCT( Kic + 2) 



(D4) 



ftic + 1 



where Ki c and a comes from the approximations 

Kf& « a Kic /cr, 



Fig. ID1I demonstrates accuracy of the approximation 
(left and middle panels show in fact the variation of I lc 
because both q and ^ ma x are constant there). Our calcu- 
lations may be summarized as follows: this approximation 
may be used, with errors less than ~ 30%, for those azimuth 
ip where e m < 1 and ef > 0.1, in the range of g from 1 — 2Ag m 
to 1, where A^ m = 1 — g m , g m is the radius (close to the 
shock) where the maximum in the radial profile of bright- 
ness happens; in addition, approximation may not be used 
for g < 0.9. If for some azimuth, the above conditions On 6m 
and ef do not hold, the accuracy of approximation gradually 
decreases because the role of the exponent in N(E) and of 
the radiative losses may not be described by the first terms 
in the decompositions used for derivation of the formula. 

Let's consider Fig. ID1I The photon energy 6m does not 
change with azimuth for IC process. On the left panels, 
the reduced fiducial energy ef(y>) = e f || /(^max^s) >0.1 for 
any azimuth: Cfy = 3 at the parallel shock and efj_ = 0.19 
at the perpendicular shock. The approximation is accurate 
for any azimuth, for 0.98<£?<lat(£ = 90° and for a 
e ff , p) (Didder range of g at if = 0°. Middle panels on Fig. ID II show 
the same case except of ef|| = 1. At parallel shock (i.e. 
ip = 90°), the range for g is smaller, 0.99 < £ < 1 (lower 
panel). Therefore, the approximation of the azimuthal pro- 
file for g = 0.98 is inaccurate (upper panel, blue line), espe- 
cially for ip > 45° where ef decreases; it is e^± = 0.06. The 
azimuthal profile is however accurate for g = 0.99 (black 
line). Similar situation is for variable i^max (right panels on 
Fig. IDl)) . ^max±/^max|| = 3.25 for considered model, there- 
fore £f||/efj_ = 52. Therefore, in order to obtain a repre- 
sentative approximation, the lowest possible Cfy should be 
about 0.1 x 52 = 5.2. We see from the figure that accuracy 
decreases toward smaller ef (i.e. where the role of radiative 
losses are very efficient in modification of the electron dis- 
tribution) and for smaller g. 

In general, the accuracy of the approximation is bet- 
ter for larger a and smaller e m . With decreasing of the as- 
pect angle <p Q , the accuracy of the approximations for the 
azimuthal profile increases at the beginning (because con- 
trasts in as, ef and £ ma x are lower) and then decreases 
again, for the case of the quasi-parallel injection, because 
SNR becomes centrally-brightened while our approximation 
is developed for regions close to the edge of SNR. 



(D3) 



© 2010 RAS, MNRAS 000.IMT51 



14 Petruk 0. et al 



p/R=0.99 
p/R=0.98 





Figure CI. Azimuthal (upper panels) and radial (lower panels) profiles of the X-ray surface brightness S x (solid lines) and its approx- 
imations (|C24|) (dahsed lines). Calculations are done for </> Q = 90°, 6 = 0, isotropic injection, 7 = 5/3, s = 2, a = 1. Models of £ max : 
£max = const (left and middle panels) and time-limited one with 77 = 1.5 (right panels). The reduced electron energy is e m = 1 and the 
reduced fiducial energy is Cfy =3 (left), efy =1 (middle), efy =5 (right panels). 




p/R p/R p/R 



Figure Dl. Azimuthal (upper panels) and radial (lower panels) profiles of the IC surface brightness S{ c (solid lines) and its approx- 
imations (|D8J) (dahsed lines). Calculations are done for cf) = 90°, 6 = 0, isotropic injection, 7 = 5/3, s = 2, a = l. Models of £ max : 
£max = const (left and middle panels) and time-limited one with 77 = 1.5 (right panels). The reduced electron energy is e m = 1 and the 
reduced fiducial energy is e-p m =3 (left), efy =1 (middle), efn =5 (right panels). 



© 2010 RAS, MNRAS 000, HHH] 



Tools for non-thermal images of SNRs 




0.92 



0.94 



0.96 



0.98 



p/R 



Figure D2. Radial profile of the radio brightness S e (solid lines) 
and its approximation (|C19J) (dashed lines) for azimuth <p = 0° 
(black lines) and (p = 90° (blue lines), <f> Q = 90°. For smaller 
aspects, <f>o < 90°, the approximation agree better with the nu- 
merical profiles. 

APPENDIX E: APPROXIMATE FORMULA 
FOR THE RADIO SURFACE BRIGHTNESS IN 
SEDOV SNR 

An analytic approximation for the azimuthal variation of the 
radio surface brightness of Sedov SNR (Paper I) may be ex- 
tended to allow also for a description of the radial variation 
close to the forward shock. Namely, the correction consists 
in a factor I r : 



where I r is the same as for the X-ray approximation (|C18[) . 
Accuracy of this approximation for the radial profile of 
brightness is demonstrated on Fig. ID2I and on Fig. 2 in Pa- 
per I for the azimuthal profiles. I r varies with azimuth less 
than 10% (cf. e.g. black and blue dashed lines on Fig. ID2|) . 
This variation is due only to /3(0 O ). Thus, f3 may be taken 
constant with a good choice 0/2 = 1 (see also Appendix |C]) . 

The smaller O , the smaller differences between the ra- 
dial profiles for azimuth ip = 0° and 90° (black and blue 
solid lines approach one another with decrease of the aspect 
angle) . 



S r (^p)0C ? (^ B M (S+1)/2 / r (p) 



(El) 
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